The radiative potential method for calculations of QED radiative 
corrections to energy levels and electromagnetic amplitudes in 

many-electron atoms 

V.V. Flambaum^ and J.S.M. Ginges^"'- 
^School of Physics, University of New South Wales, Sydney 2052, Australia 
^Department of Physics, University of Alberta, Edmonton AB T6G 2J1, Canada 

(Dated: February 9, 2008) 

Abstract 

We derive an approximate expression for a "radiative potential" which can be used to calculate 
QED strong Coulomb field radiative corrections to energies and electric dipole (El) transition am- 
plitudes in many-electron atoms with an accuracy of a few percent. The expectation value of the 
radiative potential gives radiative corrections to the energies. Radiative corrections to El ampli- 
tudes can be expressed in terms of the radiative potential and its energy derivative (the low-energy 
theorem): the relative magnitude of the radiative potential contribution is ~ a^Z^ ln(l/a^Z^), 
while the sum of other QED contributions is ~ a^(Zj + 1)^, where Zi is the ion charge; that 
is, for neutral atoms {Zi = 0) the radiative potential contribution exceeds other contributions 
~ times. The advantage of the radiative potential method is that it is very simple and can 
be easily incorporated into many-body theory approaches: relativistic Hartree-Fock, configuration 
interaction, many-body perturbation theory, etc. As an application we have calculated the ra- 
diative corrections to the energy levels and El amplitudes as well as their contributions (-0.33% 
and 0.42%, respectively) to the parity non-conserving (PNC) 6s-7s amplitude in neutral cesium 
(Z=55). Combining these results with the QED correction to the weak matrix elements (-0.41%) we 
obtain the total QED correction to the PNC 6s-7s amplitude, (-0.32 it 0.03)%. The cesium weak 
charge Qw = — 72.66(29)exp(36)theor agrees with the Standard Model value = —73.19(13), 
the difference is 0.53(48). 
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I. INTRODUCTION 



The precision of calculations and measurements of phenomena of heavy neutral atoms 
has reached the level where strong-field QED radiative corrections are observable. The most 
striking example is parity nonconservation (PNC) in the neutral cesium atom {Z = 55) 
where the nuclear Coulomb field radiative corrections "saved" the standard model of particle 
physics (this dramatic story may be found, e.g., in the review [1,]; see also the originalpapers 
for the measurement 
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While there is an abundance of highly-accurate calculations of radiative corrections to 
phenomena of single-electron or few-electron atoms, only a handful of calculations have been 
performed for atoms with many electrons. A proper account of the many-body effects in 
calculations of radiative corrections to phenomena of many-electron atoms, including the 
PNC amplitude, is lacking. 

The first estimates of radiative corrections to energies of an external electron in heavy 
neutral atoms were performed more than 20 years ago A semi-empirical formula for 
the s-wave radiative correction to energy levels (the Lamb shift) was derived. The relative 
magnitude of this correction is ~ Z'^a^\n{l/Z'^a^). It rapidly increases with the nuclear 
charge Z and was important in making a very accurate prediction of the francium {Z = 
87) spectrum (measurements performed after the theoretical prediction 16^ agree with the 
calculated energy levels to better than 0.1%). Calculations of the Lamb shift in alkali 

n 

and coinage metal atoms were performed in Ref. [17|] using local Dirac-Slater potentials; 
importantly, it was demonstrated that the ratio of energy shifts arising from the Uehling 
potential and self-energy in neutral atoms is the same as that in hydrogen-like ions, verifying 



the authors' earlier estimates 



18|. Calculations of the Lamb shift in neutral alkali atoms 



were also performed in Ref. 19^ using local atomic potentials. In our recent works ^ we 
used a parametric potential fitted to reproduce radiative corrections in hydrogen-like ions 
to perform approximate numerical calculations of radiative corrections to the energy levels, 
electric dipole (El) transition amplitudes, and a rough estimate of radiative corrections 
to the PNC amplitude in cesium. Recently, calculations of radiative corrections in local 
effective atomic potentials have been performed for El amplitudes in neutral alkalis in Ref. 
and for the PNC amplitude in Cs in Ref. |U]. 
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It is well-known that many-body effects (exchange interaction, core relaxation and po- 
larization, correlations) may be very important. For any perturbation located at small 
distances, e.g. the field (volume) isotopic shift, the matrix elements for an electron with 
orbital angular momentum / > are dominated by the many-body effects mentioned above. 
In this case one cannot guarantee the magnitude or even the sign of the shift when one 
does a model potential calculation. Moreover, even for the s-wave, many-body corrections 
can change the results by a factor of 2. QED radiative corrections also come from small 
distances, and one may expect a similar situation. Therefore, recent calculations of QED 



corrections in mode 
in a recent work 



atomic potentials cannot guarantee results of high accuracy. Indeed, 
2l| the importance of a proper account of core relaxation in calculations 
of the radiative shift due to the Uehling potential was demonstrated for neutral cesium. 

In the present work we suggest a simple radiative potential approach which allows one 
to calculate radiative corrections to many-electron atoms including many-body effects. In 
particular, our method is valid for calculations of the Lamb shift and radiative corrections 
to El amplitudes. We claim an accuracy of a few percent for calculations of radiative 
corrections to s-wave energy levels, s-p intervals, and s-p El amplitudes for neutral cesium 
using the radiative potential approach. We believe that our approach is complementary to 
the direct Feynman diagram calculations in a model potential. 

This paper is organized as follows. In Section |n] we derive an approximate ab initio for- 
mula for the radiative potential for Za <^ 1. Then we refine this potential to include higher 
orders in Za using published results for hydrogen-like ions j^]. In Section ITTTl we describe 
the procedure for calculations of QED radiative corrections to electric dipole amplitudes. In 
Section IIII Al we derive the low-energy theorem: the vertex and normalization corrections 
are expressed in terms of |^ where S is the electron self-energy operator. Further discussion 
is presented in Section IlII Bl The dominating contribution is due to the radiative corrections 
to electron wave functions produced by the radiative potential. The relative value of this 
contribution is ~ Z"^ ln{l / a"^ Z"^) both in ions and neutral atoms. It is shown in Section 
IIII (^1 that the sum of all other terms (the vertex correction, the normalization correction, 
and part of the self-energy operator which can be presented as = (if - E)A + A{H - E) 
and does not contribute to the energy shifts) is ~ <y^{Zi + 1)^ where Zi is the ion charge, H 
is the atomic electron Hamiltonian, and A is an operator defined in Section IlII CI Therefore, 
in neutral atoms {Zi = 0) this sum is times smaller than the radiative potential contri- 
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bution and may be neglected. As an application we calculate in Section IIVI the radiative 
corrections to the energy levels and El amplitudes in neutral cesium including many-body 
effects: core relaxation, core polarization by the photon electric field, and correlation correc- 
tions. In Section 13 we calculate the radiative corrections to energy levels and El amplitudes 
contributing to the parity non-conserving 6s-7s amplitude in Cs. Finally, the conclusion for 
the radiative potential approach for the calculations of the QED radiative corrections to 
energies. El amplitudes, and the PNC amplitude is given in Section IVII 

II. THE RADIATIVE POTENTIAL 

A. Derivation of the radiative potential; radiative shifts in H-like ions 

We define a radiative potential L such that its average value coincides with the radiative 
corrections to energies, 

6Er, = {n\L\n) ; (1) 

the radiative potential is non-local and energy-dependent, L = L{r,r' , E), where E is the 
electron energy. It contains the non-local electron self-energy operator in the strong Coulomb 
field, Tj{r,r',E), and the local vacuum polarization operator comprised of the lowest-order 
in Za Uehling potential $(/(r) and the higher-order Wichmann-Kroll potential. Diagrams 
for the radiative energy shifts are presented in Fig. ^ 




(a) (b) 



FIG. 1: Diagrams for radiative energy shifts, {L), corresponding to (a) the self-energy and (b) 
vacuum polarization. Double line denotes a bound electron; wavy line denotes a photon. 

The actual problem is the calculation of the self-energy S(r, r', E) contribution to the 
radiative potential. This calculation can be divided into two parts: one in which the electron 
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interaction with virtual photons of high-frequency is considered, and one in which virtual 
photons of low-frequency is considered. In the high-frequency case the external field (the 
nuclear Coulomb field) need only be included to first order (vertex diagram). In the case of 
a free electron the vertex diagram give s the electric f{q'^) and magnetic g{q'^) formfactors 
presented, for example, in the book j23|. The calculations of the contributions of /{q"^) and 
giq"^) to the radiative potential L are similar to the calculation of the Uehling contribution 
presented in js^. Therefore, we present the calculations very briefly. We also present 
the well-known results for the Uehling potential for comparison with the fiq"^) and giq"^) 
contributions calculated in this work. 

In the momentum representation the high-frequency contribution to the radiative poten- 
tial is equal to 

$rad(q) = <5rad(q)$(q) , (2) 

where $ is the atomic potential which at small distances is equal to the unscreened nuclear 
electrostatic potential, and 



Qrad(q) = - A^(-q') + 4r^7 • q + /(-q' 



(3) 



~ ' 2m 

Here the first term contains the polarization operator P(— q^) and leads to the Uehling 
potential, 7 are the Dirac matrices, and m = l/vc where Tc is the Compton wavelength. 
We use units h = c = 1 throughout, except where presented explicitly. In the coordinate 
representation 
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$rad(r) 
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where y = yq^- A method to calculate this integral is suggested in 
of P(— q^), /(— q^) and (?(— q^) from j3l we obtain 



vi.dK-y'^) ew{'iry)ydy , (4) 

2311. After substitution 



$rad(r) = $c/(r) + $g(r) + $)(r) . 



(5) 



The first term is the well-known Uehling potential 



$[/(r) = ^<l>(r) 



For the magnetic formfactor contribution we obtain 

a 



rft^^(l + -)e- 
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$(r)( / dt ; 
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(6) 



(7) 



A straightforward calculation for the electric formfactor gives 

ry /"^ lr/l\/ \S 

$J(r) = --<|.(r) )(ln(t2-l)+ln(4mVA^))-- + 



t2 



g— 2trm 



The expression for the electric formfactor © contains a low-frequency cut-off parameter 
A in the argument of the logarithm. In the standard calculation of the energy shift Q 
this parameter is assumed to be in the interval (Za)^m ^ A ^ m. After addition of the 
low-frequency contribution the parameter A cancels out. To minimize the low- frequency 
contribution we select A of the order of the electron binding energy in atoms, A ~ {Zaym, 
which is the smallest possible value that can be taken as the border for the high-frequency 
region (one can use free-electron Green's functions to calculate the electron formfactors, 
used above, for frequencies u ^ {ZaYm). This gives In(4m2/A^) = 41n(l/Za) + const, 
where the constant const does not depend on Z for Za ^ 1, const ~ 1. Therefore, the 
high-frequency contribution allows us to determine the electric formfactor contribution with 
logarithmic accuracy. The uncertainty is due to the omitted low-frequency contribution for 
this term. 

The constant in the electric formfactor can be found from a comparison with calculations 
of the Lamb shift for the high-energy levels (principal quantum number ^ 1) in hydrogen- 
like ions. Indeed, the energy of an external electron in a many-electron atom is extremely 
small, E ~ lO^^mc^. Therefore, we need the self-energy operator T,{E ^ 0). The short- 
range character of S(r, r', 0) means that we can use unscreened Coulomb Green's functions 
to calculate it. 

For very light atoms {Z ~ 1) the comparison can be made with the non-relativistic 
calculations presented in the book [23:], yielding const ^ —0.63. However, this result is not 
applicable for Z > 10. In this case we have to use the results of all-orders in Za calculations 
for hydrogen-like ions presented in Ref. [2^. To reproduce the results of {3] (and leave 
some room for a low-frequency contribution discussed below) we select this logarithm in 
the form 41n (1/Za + 0.5), where a small constant 0.5 is added into the argument of the 
logarithm. This selection gives an accuracy ~ 10% for all important applications: the s- 
wave self-energy, s-p intervals (which are needed to calculate the parity violation effects), 
and fine structure intervals for any Z. However, the radiative shifts for p-waves (and higher 
waves) are small and sensitive to the low frequency contribution. To make our calculation 
complete and improve the accuracy to 1% we should consider this contribution too. 
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A consistent calculation of the low-frequency contribution to the non-local self-energy op- 
erator S(r, r', E) using Coulomb or parametric potential Green's functions is a complicated 
task. However, at the present level of experimental accuracy this low-frequency problem is 
not of immediate importance. It is much easier to fit this small low-frequency contribution 
using a parametric potential $/(r). A typical frequency in the low- frequency contribution 
is ~ Els] therefore, the range of this potential is about the size of the Is orbital, as/Z. 
To reproduce the p-level radiative energy shifts we use the following expression for the low- 
frequency contribution 

$,(^) = _:^z^«We-^^/'^^ , (9) 

where e is the proton charge and B{Z) = 0.074 + 0.35Za is a coefficient fitted to reproduce 
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the radiative shifts for the high Coulomb p-levels calculated in 

Finally, we should introduce one more correction which becomes important for very heavy 
atoms, Z > 80. The potential (jHI) is not applicable for very small distances, r ^ Zar^. In- 
deed, we used an expression for the electric formfactor of a free electron. However, at very 
small distances the electron potential energy Za/r ^ m and we should use an expression 
for an off-mass-shell formfactor f{p,p') instead of f{—q^)- The formfactor f{p,p') leads to a 
non-local expression S(r,r',£') for r <^ ZaVc instead of the local potential (jH)). Integration 
of an electron wave function ip with a non-local operator ( / S(r, r', E)ip{r')d?r') makes the 
effective potential less singular. We take into account this fact by introducing a small dis- 
tance cut-off coefficient mr/(mr + 0.07Z^q;^). Our final expression for the electric formfactor 
contribution has the following form 

(T r°° 1r/1\/ x'^ln 

$;(r) = -A{Z, r)^<|.(r) dt -== [(l-— ) ( In {t' - l)+4 In {1/Za + 0.5))-^+-^ -^-''^ 



e 
(10) 

Here the coefficient A{Z,r) = (1.071 - 1.976x2 - 2.128x3 + 0.169x^)mr/(mr + 0.07Z^a^), 
where x = {Z — 80)a; A(Z, r) was found by fitting the radiative shifts for the high Coulomb 
s- levels calculated in (22| . 

Thus, we obtain the following expression for the complete radiative potential 

$rad(r) = <t>u{r) + <t>,{r) + ^^r) + ^r) + h'^f\r) , (11) 

where $c/(r) is the Uehling potential (jHI), ^g{r) is the magnetic formfactor contribution 
([7j), $/(r) is the high-frequency electric formfactor contribution ()10|). and ^i{r) is the low- 
frequency contribution 



To make the picture complete we added to the radiative potential a simplified form 
of the Wichmann-KroU potential (higher-orders vacuum polarization) which has accurate 
short-range and long-range asymptotics 0]: 



Accurate calculations of the Wichmann-Krol 



contribution in hydrogen-like ions for Z 



30, 40, ...100 have been performed in Ref. j2j|. To reproduce the Wichmann-Kroll s-wave 
shifts from 24 1 with a few percent accuracy one should take ^wci.^) = \^vvc^°i.'^')- The 
Wichmann-Kroll potential (fT^ gives a very small contr^ution which may be noticeable 
(~ 1%) only for Z > 80. (This confirms our conclusion [6] that the contribution of higher 
orders into the vacuum polarization potential is so small that it is unobservable in neutral 
atoms: for the cesium atom the Wichmann-Kroll potential gives a contribution to s-waves 
about 50 times smaller than the Uehling potential contribution; moreover, the Uehling 
potential itself gives only about 10% of the total radiative correction for s-waves and a very 
small contribution for higher waves.) 

Note that our fitting coefficient in the electric formfactor contribution A{Z, r) ^ 1 and 
the coefficient in the low-frequency contribution B is small, therefore the semi-empirical ra- 
diative potential ()11|) is always close to the result of the direct Feynman diagram calculation 
^ui'f') + ^gi'f') + ^f{r). It may look surprising that the radiative potential obtained in the 
approximation Za <^ 1 gives energy shifts which are close to the all-orders results. However, 
the higher-order corrections to the energy shifts are mainly due to the relativistic electron 
wave functions which we take into account exactly when calculating the matrix elements. 
Indeed, the Dirac wave function diverges at small distances r < as/Z: 

~ ^.V. ~ r-^^"^ = exp i-ZWlnr) . (13) 

The radiative energy shifts originate from very small distances r ~ = h/mc. Thus, we 
take into account the relativistic enhancement factor {ilj'^{rc)/'ip'^{aB/Z)) ~ [aB/ Zvc)^^^^ = 
(1/Z«)^^"^ = exp (Z^a^ In (1/Za)) in the Za dependence of the matrix elements when we 
use the Dirac wave functions. 

In TablelUwe compare the self-energy for 5s, 5pi/2, and 5j93/2 levels i n hy drogen-like ions 



m ny g 

f. u. 



calculated using the potential $c,(r) -|- $/(r) + $;(r) with those of Ref. j22|. It is seen that 
the radiative potential $rad reproduces the self-energy within a few percent for all Z. The 



TABLE I: Difference between the results of tlie Mohr-Kim self-energy calculations in Ref. and 
the radiative potential results ((Mohr — g — f — /)/(Mohr — g)) in % for hydrogen-like ions. 



z 


10 


20 


30 


40 


50 


60 


70 


80 


90 


100 


110 


5Si/2 


0.0 


0.4 


0.5 


0.3 


0.0 


-0.2 


-0.2 


0.0 


0.1 


0.1 


0.0 




-0.8 


-3.6 


-2.8 


-1.8 


-1.1 


-0.7 


-0.3 


0.1 


0.8 


1.8 


3.3 


5P3/2 


-2.5 


-8.3 


-8.9 


-7.3 


-5.2 


-3.1 


-1.1 


0.4 


1.4 


1.7 


0.8 



comparison is made for the highest available principal quantum number n = 5 in order to 
satisfy the condition E -C mc^ which is needed to calculate the radiative corrections in 
neutral atoms. However, in practice the results are good for any n > 1. Moreover, the 
potential $rad('") even gives the Is energy with reasonable accuracy, ~ 10%. Note that to 
calculate parity violation we mainly need to reproduce high s-level shifts in Cs {Z = 55), Tl 
{Z = 81), and Fr {Z = 87); the p-level shifts are very small and not important. 

The above ca.cu.at,o,Mand\Kose of Ref. Q) were performed in the Coulomb field 
of a point-like nucleus $(r) = Ze/r. The small correction due to finite nuclear size can 
be taken into account using integration over a realistic charge density for the nucleus, 
*&rad(r) = / "^ra"* '^'^'^'^^'^ ( | r — t' \) p{r')d^r' . The finite nuclear size contribution is suppressed 
by a small parameter r„/rc ~ 10^^. The results of our calculations for the neutral cesium 
atom presented later in this work include this correction. 



B. The radiative potential in atomic calculations 

The radiative potential we have derived from radiative shifts in hydrogen-like ions can 
be used in calculations of radiative shifts in ions and neutral atoms for all Z and for any 
number of electrons. 

Indeed, all electron wave functions with energy E <^ mc^ are proportional to the zero- 
energy Coulomb wave functions in the area r ~ r^, since the energy E may be neglected 
in the Dirac equation and the potential is unscreened in this region. Therefore, the ratio 
of the matrix elements of the radiative potential will be proportional to the ratio of the 
electron densities near the origin (at a given j,l). This is the reason why one may use 
parametric potentials fitted to reproduce Lamb-shifts in hydrogen-like ions (for principal 
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quantum numbers n ^ 1). Any potential of the range ~ Vc will give the same results. The 
radiative potentials Eqs. ()6I10|) belong to this class. This also explains the conclusion of 
Ref. that the ratio of the self-energy contribution to the Uehling contribution is the 
same in hydrogen-like ions and neutral atoms calculated in Dirac-Slater potentials. 

The magnetic formfactor potential Eq. (|7j) is a long-range one. However, it decays rapidly 
and its matrix elements are still determined by small distances r ~ aB/Z where all the wave 
functions with principal quantum numbers n ^ 1 and given j, I are proportional (since the 
energy <^ l-^'isl and may be neglected in the area r ~ as/Z). Thus, the radiative shifts 
at a given j, I are still proportional to the electron density in the vicinity of the nucleus (1/n^ 

23 show that this statement is 



in the Coulomb case). Numerical data presented in Ref. 
accurate to a few percent for n > 1 (exact for n ^ 1). 

Our semi-empirical radiative potential was derived in the field of the nucleus. In atoms 
there is an electron density contribution to the radiative potential. This can be found by 
integration of the point-charge radiative potential over the electron density (as in the finite 
nuclear size calculation). It is easy to show that this contribution is very small. The elec- 
tron density contribution is suppressed Z times relative to the nuclear charge contribution. 
Indeed, to find the energy shift we need to integrate the radiative potential with a squared 
external electron wave function ip'^lr). For the nuclear contribution ip"^ ~ ip'^lvc); for the 
electron density contribution ip"^ ~ ip'^las), since the radius of the electron charge density is 
of the order of the Bohr radius a^. For an external s-wave electron in a neutral atom the 
ratio ip'^{rc^/ip'^{aB) ~ Z. Thus, the nuclear contribution is Z times larger than the electron 
density contribution. A more elaborate estimate using the WKB semiclassical electron wave 
function and the Thomas-Fermi electron density confirms this simple estimate. 

Note that to estimate the electron density contribution to the self-energy operator we 
can use a semiclassical expression for S derived in [2^: 

2 

S(r, r', E) = In , . , ' ^, V\-e^r)) 6{r - r') . (14) 

Here V^(— e<I>(r)) = 47re^rie(r), ne(r) is the electron number density. This semiclassical 
expression is valid for r > as/Z. Again, an estimate based on Eq. (|THl shows that the 
electron density contribution is Z times smaller than the nuclear charge contribution. 

Another conclusion from Eq. (jl4p is that the self-energy operator is not sensitive to 
the energy of a valence electron in the area r ~ as/Z (the small- distance boundary of the 
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applicability of Eq. ()14p ) where e$(aB/Z) ~ {Z'^a'^)mc^ ^ E a^mc^. The logarithm in 
this area is equal to In (l/Z^a^), it is the same value that appears in the pure Coulomb case 
(c.f. Eq. (llUj) ). An estimate of the energy dependence in this area may be characterized by 
the ratio 

- (15) 



(S/E) Z2 

which is very small in neutral atoms. In ions this ratio is ~ (Zj + 1)^/Z^, where Zj is 
the ion charge (for a valence electron in an ion E ~ (Zj + lYa^mc^). We shall recall these 
conclusions during the discussion of the low-energy theorem for an electric dipole amplitude. 

C. Asymptotics of the radiative potential 

It may be useful to present long-range and short-range asymptotics of the radiative po- 
tentials. For mr ^ 1 

$.(r) ~ --T^rf In mr + 4 In Za)e-^™''$(r) . (17) 

Note that the asymptotics of the high-frequency contribution to $/(r) are presented as 
an illustration only. A correct (long-range) expression for large r is determined by the 
contribution of low frequencies. However, numerically this long-range contribution is not 
significant. Indeed, the radiative corrections to s-wave energies are proportional to \jv? {ip'^ 



near the origin) to an accuracy ~ 1% (see [22,l23)- This can be considered as an estimate 



of the contribution of the long-range tail and the energy dependence of S(r, r', E). 

The formfactor g at large distances gives a contribution which describes the interaction 



of the electron anomalous magnetic moment with the atomic electrostatic potential $ 23|: 



ah 

This long-range potential decreases faster than since the nuclear electrostatic potential 
is screened by atomic electrons. It gives an especially important contribution for orbitals 
with / > 0. The long-range character of this interaction guarantees that it is not very 
sensitive to higher order in Za corrections which are produced by the strong Coulomb field 
at r ~ r^. 
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The short-range asymptotics of the radiative potentials, mr -C 1, are the following: 

2a 

^u{r) = — ln(l/mr)<l>(r) , (19) 



(y. / \ 

$^(r) = ln(l/mr)(ln(l/mr) + ln(mVA^))$(r) , (20) 



a 

$„(r) = z — mr In (l/mr)7 ■ n$(r) . (21) 
2tt 

We see that the area mr < 1 is not important for the magnetic formfactor contribution. 
As we pointed out in the previous section the expression ()20|) for the electric formfactor 
contribution is not applicable for mr <^ Za. Indeed, the short-range asymptotics ()20|) can 
be obtained very easily using high-energy asymptotics of the vertex operator T^^{p,p',q) 
where p and p' are initial and final electron 4-momenta, and q is the photon 4-momentum. 
These asymptotics can be found, e,„.0.I„tKe case 

r'^(p,p',g)=7"exp[-^ln^(ln^ + 21n^)] . (22) 

4tt m^ m^ 

To obtain the result in the coordinate representation we should substituite q ~ 1/r. This 
gives the correction ()20|) to the Coulomb potential. However, in the area ^ p^,p''^ ^ 

r'^(p,p',g) =7^exp[-^ln^ln^] . (23) 

Ztt p' 

Here we have an explicit dependence on both electron momenta, therefore in the coordinate 
representation we obtain the non-local correction to the potential. This corresponds to the 
area where the Coulomb potential Za/r ^ m. 

For electron angular momentum j > 1/2, the short-range contributions of $t/(r) and 
$/(r) are suppressed. The dominating contribution is given by the long-range ^g{r). The 
result is also sensitive to the low- frequency contribution described by $i(r). The radiative 
shift for the pi/2 orbital is a special case. This orbital has a lower Dirac component pene- 
trating to the region r ~ r^, see Eq. fjl3p . The relative contribution of the short-distance 
area r ~ rc increases as ~ Z'^a'^. This leads to a cancellation of $/(r) and ^g{r) contribu- 
tions at Z ~ 40. For very large Z the short-range contribution dominates and the radiative 
corrections for pi/2 become comparable to that for Si/2- At Z = 110 the pi/2 shift is only 3 
times smaller than the si/2 shift. 
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III. ELECTROMAGNETIC El AMPLITUDES 



Diagrams for radiative corrections to the electric dipole (El) transition amplitude are 
presented in Fig. |21 The magnitude of different QED contributions to El amplitudes 
depends on the virtual photon gauge. Corrections to El amplitudes in hydrogen-like ions 
have been calculated with logarithmic accuracy (~ a^Z^ In (1/a^Z^)) in Ref. ^gI. It is 
pointed out there that in the Yennie gauge it is enough to take into account only those 
corrections to the non-relativistic electron wave functions produced by the non-relativistic 
radiative potential (containing S{r)). 



I 
I 



I 



I 
I 

X 



(b) 




I 
I 
I 



(d) (e) (f) 

FIG. 2: Diagrams for radiative corrections to the El amplitude: (a),(b) self-energy; (c) vertex; 
(d),(e) vacuum polarization; and (f) external photon polarization operator. Notations are the same 
as in Fig. ^ dashed line with the cross denotes the laser driven external photon field. 

In neutral atoms there is an additional small parameter suppressing the external photon 
vertex contribution (Fig. Efc)). The binding energies of the valence electron and the external 
photon frequency are extremely small in comparison with typical virtual photon frequencies. 
This makes corrections to the electric dipole operator (from the electron anomalous magnetic 
moment) and all contributions proportional to |^ negligible (see Eq. ()15p). 

Therefore, the radiative potential approach should work in neutral atoms even better 
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than in hydrogen-like ions. One has only to add the radiative potential to the atomic 
potential, calculate the electron wave functions and use them to calculate electric dipole 
matrix elements. It is convenient to perform the calculations in the length form for the 
electric dipole operator (D = er). 

In the following section f Section IIII A|) we derive the low-energy theorem for the El am- 
plitude (expressing the vertex contribution in terms of the self-energy). In Section Fill Bl we 
explain the suppression of the vertex corrections and the validity of the radiative potential 
approach. In Section Fill CI we perform a standard subtraction to remove ultraviolet diver- 
gences from the radiative corrections to the amplitude and estimate different contributions. 



A. Derivation of the low-energy theorem 

le high-frequency virtual photons the low-energy theorem follows 



For the contribution of t 
from the Ward identity 23 

r-(p,p,0) = ^, (24) 

where is the vertex operator, G^^ = 'y^p^ — m — M{p) is the inverse electron Green's 
function, and M{p) is the mass (self-energy) operator in the momentum representaion. In 
the length form for the electric dipole operator we only need the zero component of the Ward 
identity (F° = 7° — ^^^) since the electron-photon interaction in this case is described by 
-f^int ~ e'j^r^(f), where = r • E and E is the photon field. Therefore, the sum of the usual 
El amplitude and the vertex correction is 

where E = 'y'^M. Here it is assumed that the transformation of the mass operator from 
the momentum representation to the coordinate representation is accompanied by the anti- 
symmetrization of the operators r and M since, generally speaking, they do not commute 
(actually, this non-commutativity exists for the low-frequency contribution only). There is 

23| and the self-energy op- 



also some difference in definitions of the mass operator M in 
erator E used in this work (extra Dirac matrix 7*^; the matrix elements in [2^ are defined 
using ip = '?/'^7°, we use ip'^'Eip). Note that this low-energy theorem for the high-frequency 
contribution is valid in any order of perturbation theory (including all orders in Za) and 
holds for the renormalized operators (similar to the Ward identity). One should add to Eq. 
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()25|) the QED corrections to the wave functions |1 > and |2 > which are not shown there 
exphcitly. 

To prove the low-energy theorem for the low-frequency contribution we can use a non- 
relativistic expression for the QED correction to the electric dipole amplitude presented e.g. 
inRef. 

5 < l\z\2 >= / uduRefiu) , (26) 



Svrm^ JO 

, _ ^ < l\z\n >< n\pi\k >< k\pi\2 > < l\pi\n >< n\pi\k >< k\z\2 > 

+ E ,J ^ Vr^ ^Phr (28) 



n,k 



{El- En-uj){E2- Ek-Uj) 



< l\z\2 > ^ < l\pi\n >< n\pi\l > < 2\p^\n >< n\p^\2 > 

2 {Ei-En-ujf ^ {E^-E^-uf ^ > 

The typical frequency of a virtual photon is large in comparison with the difference 
of excitation energies of the valence electron. Therefore, we can simplify the the vertex 
contribution (Eq. (|2H1)) by replacing by E^^ since {E^ — En) ^ cu, and summing over n 
(closure). We repeat this procedure by replacing E^ by En and summing over k. The result 
for the vertex contribution can be presented in a symmetric form (to cancel the first order 
correction in {E^ — En)/uj and the commutator [-2,Pi]) 

I ^ < l\zpi\n >< n\pi\2 > < l\pi\n >< n\piz\2 > 

2^^ [E^-En-uf + {Ei-En-uf ^ > 

Here we neglected the small difference E2 — Ei. Now we see that the vertex contribution is 
indeed proportional io z^ (E oc '§§i ~(^~^^z|^)- Combining all terms and 

using the definition of S, 

^{Ei) = / ujdujRe}_^— , (31) 

37rm^ JO „ [El — En — uj) 

we obtain the low-energy theorem for the radiative correction to the electromagnetic ampli- 
tude: 

n hi2 — tin n ^1 — J^n 

-i<l|.i.>> (33) 

1 dH (9S 

+ -<l|fl|2>«l|-|l>H-<2|-|2». (34) 
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Note that we can extend this derivation to include the high frequency contribution. We 
just have to use the relativistic radiation operator 7^ exp (ikr) instead of the momentum 



28l 



29 



operator pi. (One should start from the relativistic expression for the amplitude 
instead of Eqs. ^-^.) 

In this derivation we used the long-range character of the electric dipole operator r. 
Indeed, we assumed that the matrix element < n\z\k > in the sum over n and k is dominated 
by states with \En — Ek\ ^ cj. This is certainly correct for the states located at distances 
comparable to the radius of the valence electron where \En\ ~ \Ek\ ~ \Ei\ ~ \E2\. For 
short-range operators (e.g. weak or hyperfine interactions) the contribution of states ~ 
\Ek\ ^ \Ei\ ~ |i?2| may be important. 



B. Enhancement of the self-energy contribution 

The contribution of the electron self-energy E to the El amplitude (Eq. ()H2|l ) is enhanced 
by the small energy denominator Ei^2 — En corresponding to the excitation of an external 
electron. The vertex (jH^ and normalization (jMj) contributions are not enhanced since 
1^ ~ T./uj, where to ^ £'i_2 — -E„ is a typical virtual photon frequency. One may conclude 
that the vertex and normalization contributions are suppressed relative to the self-energy 
contribution by a small factor (-Ei,2 — En) /to. Moreover, the vertex and normalization 
contributions usually have opposite signs and partially cancel each other. This is seen if 
we introduce complete sums X^n I"- >< ''^1 between the operators D and |^ in the first and 
second terms in Eq. ()33|) : the large diagonal contributions in these sums {\n >= |1 > in the 
first term and \n >= |2 > in the second) cancel exactly the normalization contribution Eq. 

(El- 
There is another way to explain the suppression of the vertex contribution. The product 
^ as(r,r^,E) gj2iall everywhere inside a neutral atom. The matrix element of r typically comes 
from the distance r ^ as- At this distance the nuclear Coulomb field is screened, Z^ff ~ 1; 
therefore, in this region ^^^g^ '^^ cannot have enhancement. On the other hand, at small 
distances where the nuclear charge is unscreened, r is small. This can also be seen from the 
vertex diagram itself where the operator r is locked inside the virtual photon loop located at 
a small distance from the nucleus. There is no such suppression for the radiative correction 
to the electron wave function. The radiative potential changes the energy of the electron. 
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This in turn changes the large distance asymptotics of the electron wave function and the 
matrix element of r. 

In the approximate expression for the radiative potential we use in this work 9nry,E) ^ Q 
and so there are no vertex or normalization contributions. The first two terms in the radiative 
correction Eq. ()32|1 can be presented in terms of corrections to the wave functions produced 
by S: 

^ <l\j:iE,)\n><n\ ^ <2\J:{E,)\n><n\ 
|is>-2^ p p ' Ms >- 2^ p p , I35j 

< 1|D|2 >s=< ls|£'|2 > + < l|£'|2s > . (36) 

We have checked that our approximate expression for the radiative potential gives correct 
diagonal matrix elements of S (radiative shifts). Here we need the non-diagonal matrix 
elements < l,2|S|n >. However, the main contribution to the matrix element of = 
er ~ eas is given by the low-energy states with En ~ ti^/{ma%) ~ £'1^2 ^, where to is 
the virtual photon energy. Therefore, our approximate radiative potential should give such 
matrix elements < 1, 2|E|n > correctly. Note, once again, that this statement is incorrect for 
the radiative corrections to the matrix elements of short-range operators like the weak and 
hyperfine interactions. In this case the states with En ~ he/ rtypicai > oj give a significant 
contribution and use of the radiative potential designed to give matrix elements between 
low-energy electron states is not justified. 

The derivation (and result) of the low-energy theorem we have presented above is similar 

to that of the low-energy theorem for correlation corrections to the electric dipole amplitude 

i 

in our work [31]. The vertex (structural radiation) and normalization correlation corrections 
are also proportional to |^ and suppressed by a factor i?vaience/-Ecoro ~ 1/10 where i^vaience 
and Ecore are ionization energies for the valence and core electrons. It is interesting to note 
that for the correlation corrections the vertex and normalization contributions were found to 
be numerically small for both long-range (electric dipole) and short-range (weak, hyperfine) 
operators. 

There is also a certain similarity between this theorem and the low-energy theorem (the 
Low theorem) for bremsstrahlung (see, e.g., [23 1). The main radiation comes from the 
external particle ends of the scattering diagram and is expressed in terms of the elastic 
scattering amplitude. The structural radiation (from inside the scattering vertex) is small 
and is expressed in terms of the derivative of the elastic amplitude. In our case we consider 
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the radiation of a weakly bound electron {E < 0) which is not so different from the radiation 
of an unbound particle {E > 0) if the energy E is small. 



C. Estimates of different QED corrections 

All terms in Eqs. ()26I28I27I29|) and ()32l33l34p are ultraviolet divergent as k tends to 
infinity. Therefore, we have to perform a subtraction of a standard counter term in the 
expression for the self-energy operator, 

Pi\n >< n\pi _ pi\n >< n\pi _ pi\n >< n\pi{E - E^) , . 

{E-En-uj) {-uj) ~ uj{E-En-uj) ' ^ ^ 

to cancel the linear divergence in Eqs. (j26|27|) and regroup other terms to cancel the loga- 
rithmic divergences. After the subtraction and commutation Pi{E2 — H) = —[pi,H] + {E2 — 
H)pi = 'iVj(— e$) + {E2 — H)pi the first term in Eq. (jTTj) can be transformed as follows (we 
use the operator form of Eq. ()27|) for brevity) : 

+ i-(.p,-.|2><2|p.)^^-]^p,. (39) 

The second term in Eq. ipTj) gives a similar contribution. 

All terms containing zVi(— e$) are combined to give a low-frequency contribution to the 
radiative potential. In particular, the large-cj contribution in Eq. ()38p gives a well-known 
local term in S 

S(r, r', E) = In V'{~e^r)) 6{r - r') (40) 

which cancels the low-energy cut-off parameter k in the high-frequency contribution propor- 
tional to In ^ (from the formfactor /; see |22|). Here V^(— e<l'(r)) = AiTe'^Z6{r). It is easy 
to estimate the contribution of S to the QED radiative corrections to < s\z\p > (see Eq. 
()32|)). The non-relativistic density of an s-wave valence electron is ^'^(0) ~ {Zi + l)^Z/a^, 
where Zi is the ion charge {Zi = for a neutral atom and = Z — 1 for a hydrogen-like 
ion), the energy Ei ~ {Zi + l)^a;^mc^. Therefore, the relative value of the QED correction 
produced by S is ~< > /Ei ~ a^Z^ In (l/a^Z^). 

The terms in Eq. fl39|l are combined with the vertex and normalization contributions to 
produce equations which do not contain ultraviolet divergences. The terms proportional to 
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zpi (see Eq. fl39p ) and piZ combined with the vertex contribution ()28|) give the following 
QED radiative correction to < l\z\2 >: 

1 1 cx X — > /I 11,11,1/,! 11 En — El -\- El, — E-) , , Er, — El , , . 

The energy-dependent factor in this equation is always of the order of unity (it varies from 
2 to Inl-f^l). The ratio vie for a valence electron is p/mc ~ iZi + l)a. Therefore, the 
relative correction from the vertex term is ~ (Zj + A more sophisticated estimate 

based on closure gives the same result: a < l\pzp\2 > /m? = a < l\p'^z + zp'^\2 > /m? = 
2a < l\{Ei + e<^)z + z{E2 + e<l')|2 > /m ~ a'^{Z, + l)/m ~ a^{Zi + l)^ < l\z\2 >. For 
hydrogen- like ions this correction is Z^a^, comparable to the radiative potential contribution 
(~ a^Z'^ In {l/a^Z'^)). However, for neutral atoms this correction is ~ a^, i.e. it is extremely 
small. 

The terms proportional to z\2 ><2\pi and Pi\l >< l\z combined with the normalization 
contribution ()29|) give the following QED radiative correction to < l\z\2 >: 

< 1\Z\2 >,orm= -^;^{l\z\2)[{l\p^\l) + (2|/|2)] . (42) 

For a valence electron p/m ^ {Zi + l)ac, therefore the relative QED correction can be 
estimated as ~ {Zi + l)^a;^, of the same order as the vertex contribution. 

Finally, the contribution of the external photon polarization operator (Fig. I2tf)) also 
does not have the enhancement, since it comes from r ~ a^. 

Thus, the only important contribution in neutral atoms is that of the radiative potential. 



IV. APPLICATIONS TO NEUTRAL CESIUM 



In this section we apply the radiative potential method to calculate radiative corrections 
to energy levels and electromagnetic amplitudes in the neutral cesium atom. We limit our 
consideration to the s and p levels of the external electron which are important for the parity 
violation calculation. All calculations are performed taking into account finite nuclear size. 



A. Energies 



To calculate the radiative corrections to the energy levels we add the radiative poten- 
tial to the nuclear Coulomb potential and calculate the self-consistent direct and exchange 
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potentials obtained using Dirac-Hartree-Fock (DHF) equations for the electron core. This 
DHF potential includes the potential 6V ("core relaxation") which arises from the change 
in the core electron wave functions due to the radiative potential. Then we calculate the 
energy levels of the external electron in this DHF potential produced by the core electrons. 



The next step is to include the correlation corrections. 
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t is convenient to calculate these 
33l which takes into account all 



corrections using the correlation potential method 
second-order correlation corrections and three dominating series of higher-order diagrams 
(screening of the electron-electron interaction, the hole-particle interaction, and iteration of 
the correlation self-energy) to all orders in the residual Coulomb interaction. The non-local 
and energy-dependent correlation potential E™'"^(r, r', E) is defined by the equation for the 
correlation correction to the electron energy 6E!^°^^ = {n\T,'^°^^{r, r', E)\n); it is defined in an 
analogous way to the radiative potential, Eq. (0). We add to the Dirac-Hartree-Fock 
potential to include it to all orders. The results of calculations for the radiative corrections 
are presented in Table ITTl 



TABLE II: Radiative corrections to ionization energies in Dirac-Hartree-Fock field without relax- 
ation ((DHF)o), with relaxation included ((DHF)o + (5^), and with correlation corrections included 



((DHF)o + 6V + S^"'-^); units cm 



-1 



level 


6Sl/2 


7si/2 


6pi/2 


7pi/2 


8pi/2 


9pi/2 


(DHF)o 


15.5 


4.3 


0.2 


0.07 


0.03 


0.02 


(DHF)o + dV 


15.9 


4.3 


-0.8 


-0.3 


-0.1 


-0.07 


(DHF)o + SV + E™'^'' 


17.6 


4.1 


-0.4 


-0.1 


-0.05 


-0.03 



It is seen that the many-body corrections change the result for the s-levels by ~ 10% and 
they change the sign and magnitude for p-levels. Our results for Uehling relaxation (not 
explicitly presented) are in perfect agreement with those of Ref. 2l| . 

The results are in agreement with our previous calculations jC | and lie within the range 



spanned by model potential calculations of the 6s Lamb shift performed in Ref. 



17[ (from 



15 to 27 cm ^) and Ref. (from 13 to 23 cm ^). However, the accuracy of our present 
work is higher, ~ 1% for s-levels. 
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B. El amplitudes 



We use a similar method as that for the energies to calculate the radiative corrections to 
the electromagnetic amplitudes between the s and p levels. First, we calculate the external 
electron wave functions including the radiative potential and the core relaxation 6V. Then 
we use these wave functions to calculate the radiative corrections to the electromagnetic 
amplitudes in the Dirac-Hartree-Fock approximation. At the second step we calculate the 
effect of the electron core polarization by the photon electric field using the time-dependent 
Hartree-Fock method. These core polarization corrections are often called the RPAE (ran- 
dom phase approximation with exchange) corrections. At the final step we use the correlation 
potential method j^, 33| to calculate the correlation corrections to the radiative corrections. 



In fact, to the required accuracy it is enougth to add S^'"''" to the Dirac-Hartree-Fock equa- 
tions and calculate the external electron wave functions. Other correlation corrections are 
proportional to and contribute about 1% only (the low-energy theorem; we mentioned 

this at the end of Section QllBD . The results of the calculations for the radiative corrections 

n 

are presented in Table lllll Following Ref . j20| , we present the results in terms of the di- 
mensionless relative radiative corrections R to the electromagnetic amplitude defined by the 
relation 

(y 

{s\r\p) = {s\r\p)o{l + -Rsp) . (43) 



TABLE III: Relative radiative corrections Rsp to the electromagnetic amplitudes (si/2k|Pi/2) = 
(■5i/2k|Pi/2)o(l + ^^sp) in Dirac-Hartree-Fock field with relaxation included (DBF), 
with RPAE corrections included (DHF+RPAE), and with correlation corrections included 
(DHF+RPAE+S'=°''''). 



Transition 


6s — 6p 


7p 


8p 


9p 


7s — 6p 


7p 


8p 


9p 


DBF 


0.266 


-2.90 


-4.62 


-5.68 


-0.451 


0.270 


-2.07 


-3.21 


DBF+RPAE 


0.286 


-4.39 


-11.9 


-29.7 


-0.432 


0.270 


-2.20 


-3.60 


DBF+RPAE+S^°"' 


0.265 


-2.91 


-6.25 


-10.3 


-0.340 


0.231 


-1.60 


-2.52 



Note that the radiative corrections for the 6si/2 — 7pi/2, 8pi/2, 9pi/2 amplitudes are large 
since these amplitudes are small and sensitive to any corrections to the DHF potential. 
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Unfortunately, we calculated but did not present the radiative corrections to the electro- 



magnetic amplitudes in our previous paper 



we presented only their total contribution to 



the parity violating amplitude. Anyway, the accuracy in our present work is higher (~ 1%). 



In a recent work 20|, direct calculations of the radiative corrections to electromagnetic 
amplitudes in the neutral alkali atoms were performed. In particular, the radiative correction 
to the 6s — 6p amplitude in Cs was calculated using the local Kohn-Sham potential. In this 
approach, the non-local exchange interaction is replaced by a semi-empirical local term that 
depends on the electron density. Ref. j^l also does not take into account many-body 
corrections: core relaxation, RPAE, and correlation corrections. Nevertheless, 6,s — 6p is 
a large "resonance" amplitude, and the many-body effects should not be very significant 
here. Indeed, our Dirac-Hartree-Fock value R = 0.266 is very close to the result of Ref. |2oi |: 
R = 0.261. We select our Dirac-Hartree-Fock value for comparison since it is the lowest- 
order approximation we use (no RPAE or correlation corrections) and most similar to that 
of Ref. |20| . However, the very small difference between the results is accidental. The result 
of Ref. 120] does not include the Uehling potential contribution (~ —10%), and one may 
also expect some difference due to the core relaxation effect and the different treatment of 
the exchange interaction. 

V. RADIATIVE CORRECTIONS TO THE PNC AMPLITUDE IN CESIUM 

Now we can calculate the QED radiative corrections to the PNC 6s — 7s amp litude in 
Cs. It is convenient to use the results of the sum-over-states approach of Ref. jSJ]: 

{7s\D\np) {np\Hw\6s) {7s\Hw\np){np\D\6s) 
EpNc = 2^ p _ p + p _ p 

= (1.908 - 1.352 - 0.070 - 0.020 + ...) + (-1.493 + 0.120 + 0.010 + 0.003 + ...) 
= -0.894 + ... . (44) 

The units are ieaB^—Qw/N) x 10^^^, where Qw is the nuclear weak charge and is the 
number of neutrons. 98% of the sum is given by the terms n = 6, 7, 8, 9 explicitly presented 
above. The result includes the many-body corrections to all matrix elements. 
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A. Contributions of energies and El amplitudes 



Now it is very easy to calculate the contributions to the PNC amplitude from radiative 
corrections to the energy intervals and electromagnetic amplitudes. Using the last line of 
Table IH] (with all many-body corrections included) we find that the radiative corrections to 
the energy intervals change the PNC amplitude by -0.33%. Using the last line of Table ITTTl 
(with all many-body corrections included) we find that the radiative corrections to the El 
amplitudes {Qs,7 s\D\npii2) change the PNC amplitude by +0.42%. 

Note that in Cs these two corrections nearly cancel each other, the sum of the two 
contributions is A=(0.42-0.33)%=0.09%. For the first time this cancellation was noted in 



our work 



where these corrections in the DHF approximation were estimated to be 0.33% 
(El) and -0.29% (energies). We found that the difference between our old and new results is 

n 

mainly because in p we used only 3 dominating terms in the sum ()44|) while in the present 
work we use 8 terms. The fact that in we used a different (and less accurate) radiative 
potential is not so significant. The contribution of the radiative corrections to the omitted 
terms n > 9 in the sum (j44j) in the present work is estimated (using an asymptotic formula) 
to be ~ 0.01%. 

To estimate the error we also performed another calculation. We neglected the low- 
frequency contribution $i(r) and set the coefficient A{Z,r) = 1 in $/(r). This variation 
changes the s-wave shifts by a few per cent only. However, the p-wave shifts change several 
times since they are small and sensitive to the low- frequency contribution $/(r). In this 
case the El contribution is 0.44%, the energy contribution is -0.35%. However, the sum A= 
0.09% does not change. Thus, the value of A is very stable and practically does not depend 
on the choice of the (short-range) radiative potential if this potential gives correct energy 
shifts. We estimate the uncertanty in A = 0.09% as ~ 0.01%. 

B. Weak matrix elements 

The self-energy and vertex QED radiative corrections to the weak matrix elements have 
been calculated in Refs. [s , 0, [l^ Illl, IlE 3) using Coulomb wave functions. However, all 
neutral atom wave functions near the nucleus are proportional to the Coulomb wave func- 
tions since the screening of the nuclear Coulomb potential and the electron energy may be 
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neglected here. Therefore, the relative magnitude of the QED radiative corrections for an 
external electron in a neutral atom is the same for all weak matrix elements {npij2\Hw\n' s) 
and coincides with the Coulomb case for n, n' ^ 1. Moreover, even the many-body correc- 
tions do not influence this statement since these corrections are proportional to the weak 
matrix elements {npii2\H\Y\n' s) and are multiplied by the same factor (equal to the rel- 
ative QED correction for the weak matrix element) which does not depend on n and n' . 



B,y,0, n, 12, Q 



to find the 



(lowest order); -0.9(1) |ll| (low- 
Il2j (lowest order and higher orders in 



This means that we can use the Coulomb results of Refs 
contribution of the QED corrections to the weak matrix elements. The results of calcu 
lations for the self-energy and vertex contribution are the following (in %): -0.73(20) 
(all orders in Za, using approximate relation); -0.6 
est order and estimate of higher orders )^0. 85 
the logarithmic approximation); -0.815 13| (all orders 2s — 2p). Note that even for n = 2 
the binding energy to mass ratio is En/mc^ = 0.02 and the relative radiative correction 
should be very close to the large n limit (i.e. one may expect a difference with the large 
n limit result ~ En/mc^ = 2%). Based on these results we assume the correction equal to 
(—0.82 ± 0.03)% which is in agreement with all calculations. 



The Ue: 
to Refs. 



me 



)otential contribution is easy to find and the error is negligible. According 
the Uehling contribution to the weak matrix element is 0.42% 



(note that the Uehling contribution to the weak matrix element is practically the same as 
its contribution to Ep^c since the sum of the Uehling contributions from the energies and 
electromagnetic amplitudes is about -0.01%). The Wichmann-KroU contribution is very 
small, -0.005% QQ. 

Thus, the sum of all radiative corrections to the weak matrix elements contributes 
(-0.82 + 0.42 - 0.005)% ^ (-0.41 ± 0.03)% to Epnc- 



C. PNC amplitude 

The sum of all QED radiative corrections to Ep^c is -0.33% (energies) -|- 0.42% (El) - 
0.41% (weak) = -0.32 ± 0.03%. The error ~ 0.01% in A =-0.33% (energies) + 0.42% (El) 
gives a small contribution to the error in Ep^c if added in quadruture. Note that all three 
contributions (energies. El, weak) are equally important here, and the El contribution is 
the largest. 
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Recently, calculations of radiative corrections to the PNC amplitude in Cs were performed 



in Ref. in an effective atomic potential. The result is (—0.27 ± 0.03)%. Our result 
(— 0.32±0.03)% is slightly different possibly due to the many-body corrections which have not 
been calculated in the work 14] . Note that the net effect of the many-body corrections in Cs 
is not very large because of the accidental cancellations of different many-body corrections. 
One can see these cancellations between the RPAE and correlation contributions in Table 
for the radiative corrections to the El amplitudes. A strong accidental cancellation happens 
also for the main PNC amplitude where there are 4 large corr_e 
each, and the sum of all correlation corrections is 2% only 
not take place, for example, in the PNC amplitude for Tl, where the total contribution of 
the many-body corrections is large. 

The many-body calculations of the PNC amplitude produced by the electron-nucleus 



ation corrections, up to ±20% 
, 0| ; such a cancellation does 



weak interaction are described in detail in our review 



where one can a 



references (see also the original papers for calculations of atomic structure 



so find numerous 

^ . and 

Breit corrections ^5lh6lh3,h8|). The contribution of the weak electron-electron interaction 

fin 

is very small. Within the standard model it is equal to 0.04% p, Therefore, the PNC 
amplitude is proportional to the nuclear weak charge Qw The result, including the radiative 
correction (—0.32 ± 0.03)% calculated in the present work, is the following: 

EpNc = -0.898(1 ± 0.5%) X 10-^HeaB{-Qw/N) . (45) 



From the measurements of the PNC amplitude we obtain 

Qw = -72.66(29)exp(36)theor • (46) 
The difference with the standard model value = —73.19(13) is 

Qw-Q^w= 0.53(48) , (47) 

adding the errors in quadrature. 

VI. SUMMARY AND CONCLUSION 

We suggest to calculate radiative corrections to energy levels and electromagnetic ampli- 
tudes using the simple radiative potential 

$rad(r) = $c/(r) + $,(r) + $^(r) + $,(r) + ^$rc"(^) , (48) 
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where ^ui'f') is the Uehhng potential ^gi'f') is the magnetic formfactor contribution 
$/(r) is the electric formfactor contribution (fTUI) . $/(r) is the low-frequency contribution 
(jni). The simplified Wichmann-KroU potential $^^'*^(r) (jl2j) gives a very small contribution 
which may be noticeable (~ 1%) only for Z > 80. 

The results obtained using this radiative potential are in good agreement (few percent) 
with the radiative corrections to the s, pi/2 and P3/2 energy levels calculated in [2^ for the 
hydrogen-like ions. The advantage of the radiative potential method is that it is very simple 
and can be used in many-electron atoms and molecules. 

To calculate radiative corrections to electric dipole amplitudes we suggest use of the 
low-energy theorem derived in Section Fill Al It is applicable because the ionization energy 
of an external (valence) electron E^x and the external photon frequency cUex are small in 
comparison with the typical frequency of a virtual photon uj^; Eex, t^ex ^ The vertex 
and normalization corrections are expressed in terms of the energy derivative of the electron 
self-energy operator. The dominating contribution is given by the corrections to the elec- 
tron wave functions produced by the radiative potential. The relative contributions of the 
remaining corrections are small, ~ in neutral atoms and ~ {Zi + 1)^/Z^ in ions. 

The radiative potential method allows us to take into account many-body effects. First, 
we add •I'rad to the Dirac-Hartree-Fock Hamiltonian (i.e., use the potential V = Vdhf + ^i:a,d) 
and calculate the new self-consistent field which includes a correction to the Dirac-Hartree- 
Fock potential SVdhf due to the change of the internal electron orbitals produced by $rad 
(the relaxation effect). The relaxation effect is always significant. For a short-range potential 
$rad it is larger than the direct radiative shift for electron angular momenta j > 1/2. Then 
we can use the new electron energy levels and DHF orbitals to calculate the correlation 
corrections applying the many-body theory methods which are described, for example, in 
Ref. y. 

We applied the radiative potential method to calculate the radiative corrections to en- 
ergy levels and electromagnetic amplitudes in the neutral Cs atom and demonstrated the 
importance of many-body effects in such calculations. The many-body effects change the 
s-level radiative shifts by ~ 10%, and they change the sign and magnitude of the p-level 
shifts. Many-body effects in the radiative corrections to the electromagnetic amplitudes 
are also very significant. The RPAE (core polarization) corrections usually enhance the 
radiative correction. The effect is especially significant for the small amplitudes where we 
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observed the RPAE enhancement up to 6 times. The correlation corrections usually act in 
the opposite direction and are equally significant. 

Finally, we calculated the contributions of the radiative energy shifts and radiative cor- 
rections for the electromagnetic amplitudes to the PNC 6s — 7s amplitude in cesium. The 
radiative corrections to the weak matrix elements have been calculated in previous works. 
The sum of all QED radiative corrections to Ep^c is -0.35% (energies) + 0.42% (El) - 0.41% 
(weak) = (— 0.32±0.03)%. Note that all three contributions are equally important here, and 
the El contribution is the largest. Using this radiative correction and previous many-body 
calculations we obtain the PNC amplitude Ep^c = — 0.898(1±0.5%) x lO'^^ieapi—Qw/N). 
From the measurements of the PNC amplitude 0] we extract the Cs weak charge = 
— 72.66(29)cxp(36)theor- The difference with the standard model value = —73.19(13) 
Q is Qh^-QSM^ 0.53(48). 
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